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Abstract 

The  earth  is  usually  represented  in  Cartesian  Coordinates  on  a  projection 
plane  with  the  atmosphere  mapped  three-dimensionally.  Most  such  maps  are 
orthogonal  and  in  this  case  the  equations  of  motion  take  a  particularly  simple 
form. 
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Transformation  of  the  Equations  of  Motion  of  Meteorology 
into  Arbitrary  Orthogonal  Coordinates 


1.  INTRODUCnON 

The  flow  at  various  levels  in  the  atmosphere  is  generally  represented  on 
flat  surfaces,  or  maps.  The  equations  which  govern  this  flow  are  derived  for 
the  spherical  earth.  In  order  to  solve  these  equations  numerically,  as  applied 
to  actual  maps,  one  must  transform  the  flow  equations  into  map  coordinates. 

The  coordinate  system  to  be  used  depends  upon  the  problem.  For  example, 
almost  all  numerical  prediction  work  concerns  itself  with  flow  over  a  large  part 
of  the  Northern  Hemisphere  represented  on  a  stereographic  projection.  For 
some  purposes,  the  Mercator  projection,  or  a  combination  of  Mercator  and 
stereographic,  have  been  used. 

For  easy  reference,  it  seemed  desirable  to  write  the  equations  of  motion  in 
several  orthogonal  coordinate  systems.  In  this  paper  the  equations  of  motion  of 
meteorology  have  been  transformed  from  spherical  coordinates  into  any  orthogo¬ 
nal  system,  so  that  one  need  only  substitute  the  transformation  formulas  without 
actually  carrying  out  the  lengthy  computations  to  obtain  the  desired  system. 
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2.  THE  TRANSFORMATIONS 


The  equations  of  motion  of  meteorology  are  usually  expressed  in  terms  of 
distances  and  directions  relative  to  the  earth.  Thus,  we  have: 


Fg  ■  r  cos  0X  +  2  (X  +  fl )  (  r  cos  0  — r  sin  0  <t> ), 

Fft  ■  r  it  2r  ^  +  r  sin  0  cos  0  X(X  +  2  0),  (1) 

Fu  *  r  — r  0^ — r  cos^  0X  (X  +20). 


Here,  Fg,  Fr,  Fu  are  respectively  the  eastward,  northward,  and  upward  com¬ 
ponents  of  the  force  on  a  particle  of  unit  mass;  X.,  0,  r  are  the  usual  spherical 
coordinates— longitude,  latitude,  and  distance  to  the  center  of  the  earth,  respec¬ 
tively— and  O  is  the  angular  velocity  of  the  earth's  rotation;  the  dots  represent 
differentiation  with  respect  to  time. 

However,  in  many  situations  in  which  the  equations  are  to  be  applied,  one  is  con¬ 
cerned  not  with  the  spherical  earth  but  with  a  flat  projection  or  map  of  the  earth's 
surface,  the  atmosphere  being  represented  three-dimensionally  above  the  projection. 
For  such  applications  it  may  be  useful  to  have  an  equivalent  set  of  equations  which 
involve  only  quantities  in  the  projected  system.  The  purpose  of  this  report  is  to 
demonstrate  how,  for  a  large  class  of  maps,  the  equations  of  motion  can  be  trans¬ 
formed  into  a  simple  equivalent  set  of  equations  relative  to  the  map. 

We  consider  the  earth  as  a  sphere  of  radius  r0>  X.,  0,  r  are  spherical  co¬ 
ordinates  as  above.  Then,  if  the  projection  plane  i6  the  x-y  plane  of  a  rectangu¬ 
lar  system  x-y-z,  the  map  is  determined  by  expressing  x  and  y  as  functions  of 
X  and  0: 


x  *  x  (X,  0), 
y  *  y  (X,  0) 


(2) 


and  z  as  a  function  of  r  alone,  specifically: 


o* 


z  ■  r—  r 


(3) 
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chosen  so  that  z  •  o,  that  is  the  x-y  plane,  corresponds  to  the  earth's  surface 
r  *  r  .  The  inverse  transformation  is  then  representable  in  the  form: 

i  *  4  (x,  y), 

X  «  X  (x,  y),  (4) 

r  *  2  +  r  . 
o 

AAA 

Now  let  i,  j,  k  be  the  unit  vectors  in  the  spherical  system  with  directions: 
l  -  eastward, 

A 

j  -  northward, 

A 

k  -  upward. 

Velocity  and  force  vectors  may  then  be  expressed  in  terms  of  these.  The 
velocity  V  is  Been  to  be: 

V  »  r  cos  $xt  +  r$J  +  r  k,  (5) 

and  the  force  F  on  a  unit  mass,  is  simply: 

F  *  F  f  t  F  ]  +  F  5, 

eA  nJ  u  ’ 

(6) 

Fe'  Fn'  Fu  as  111  E<1* 

The  position  vector  R  (relative  to  the  center  of  the  earth)  is: 

?T«  r  k  (?) 

The  partial  derivatives  of  this  expression  with  respect  to  x,  y,  z  are  vectors  in¬ 
dicating  the  directions  of  increasing  x,  y,  z  respectively: 


-L<r£).r8l.  Yafcax*  + 

8x  'rK'  r  8  x  r\8X  8x  8  ^  8  x  +  8r  8  x  / 


and  since 


„„  ||.o. 
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In  the  following  we  shall  assume  that  this  orthogonality  condition  holds  in  the 
given  projection. 

It  is  now  possible  to  derive  expressions  for  the  components  of  the  force  and 

the  velocity  in  the  new  system  by  dot  multiplying  the  respective  vectors  with  the 

unit  vectors  R  ,  R  ,  R  .  Thus,  let  F  .  F„,  F_  be  the  components  of  the  force  in 
x  y  z  x  y  z 

the  x,  y,  z  directions  respectively,  and  let  u,  v,  w  be  the  components  of  the 
velocity. 
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F  cos  a  +  p  JtA. 
-•a  *e  9  3x  n  3x 

fx  .  F  •  rx - r - 

.  F  cos  ^  4“  +  F  4^ 
F  .?.£ 


Fz  ■  F  •  Rz  *  Fu- 


u  “  V1-  Rx  =  t ^  (cos2  +|4*)' 

*  v"' VtH2  *t?X  +lf  0' 


w  ■  v“.Aa-r. 


Solving  for  X  and  $  in  the  first  two  equations  of  Eq.  (14)  we  get: 


ut  ii_  vt  aj 


r  cos  4  J 
ut  3  X  vt  31 

xyy-  yy? 


where  J  is  the  Jacobian  of  the  transformation  (X,  $  — ►(x,  y),  that  is. 


ax  a»  a x  m 

J  SxTy  Ty  3  x  1 


Substituting  Eq.  (15)  in  Eq.  (13)  and  using  Eq.  (1),  Eq.  (14)  and  the  condition 
Eq.  (12)  one  obtains: 


3 


_  du 
F  ■  X-  —  v 
x  at 


v&*“ 


4 


dw  (u2+  v2) 


[—  *-^v  (■£-'&)]” 

-z  "ii- lir^J-20cos2  * 

where  e  »  sign  J 

(*■“*■  •■!:!!!}?» 

Equations  (17)  are  the  equations  of  motion  in  the  projected  system. 


(17) 


(18) 


3.  NOTE  ON  THE  ORTHOGONALITY  CONDITION,  EQ.  (12) 


From  a  brief  survey  of  the  maps,  or  projections,  in  common  use,  it  is  seen 
that  this  condition  is  a  reasonable  restriction  as  for  most  such  maps  it  is  already 
satisfied.  As  a  matter  of  fact,  one  of  the  most  common  restrictions  map-makers 
wish  to  impose  is  the  more  general  one  of  conformality;  that  is,  preservation 
of  angles  or,  equivalently,  equality  of  scale  in  all  directions  at  a  given  point 
(actually,  then.the  preservation  of  shapes  of  small  regions).  Analytically  the 
conformality  conditions  may  be  expressed  by  Eqs.  (19a)  and  (19b): 


8k  8k 
8x  By 


(t  t  as  above). 

x,  y 


(19a) 

(19b) 


Eq.  (19a)  is  of  course  the  orthogonality  condition  Eq.  (12),  and  so  conformality 
implies  orthogonality  (but  not  conversely  because  of  the  additional  condition 
Eq.  119b]  ) . 

One  further  remark  on  condition  Eq.  (12).  Itcanbe  expressed  equivalently  by 
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Indeed,  it  ie  perhaps  more  useful  in  this  form,  since  one  is  usually  given  x 
and  y  in  terms  of  X  and  4  rather  than  the  inverse  direction  which  would  have  to 
be  computed  before  condition  Eq.  (12)  could  be  applied. 


4.  THE  TRANSFORMED  EQUATIONS  FOR  SOME  SPECIFIC  MAPS 


One  of  the  maps  most  in  use  today  is  the  Stereographic  Projection.  This  con  - 
sists  of  a  projection  of  the  earth's  surface  upon  some  secant  plane  from  the  point 
on  the  earth  opposite  the  plane.  If  we  take,  as  is  usually  the  case,  the  projection 
plane  perpendicular  to  the  earth's  axis  and  the  point  of  projection  as  the  south 
pole,  then  by  geometry  we  can  compute  the  equations  of  the  map. 


Then  we  have 


c  COB 

(1+  sin 


£> 


y 


C  COS  4 

(1  +  sin4 ) 


cos  X, 
sin  \, 


(21) 


z.  ■  r  -  r 
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and  the  inverse  transformation 


X  «*  tan  *(y/x)  , 


sV 

\c2  +  <x2  +  y2)/ 


(22) 


r  -  z  +  rQ  . 


The  orthogonality  condition  is  now  easily  verified  in  Eq.  (12)  or  Eq.  (20). 
(Actually,  this  is  unnecessary  since  it  is  well  known  that  the  stereographic  projection 
is  conformal  and  therefore  orthogonal.) 

Applying  these  equations  to  Eq.  (17)  we  get  for  the  stereographic  projection; 


F  -  — - v 
x  dt 


Fy  •■ai+  u 


/  2  /  2  2  a  n 

/c  -(x  +  y  )\  1 

(uy— vx) 

4-  ur 

4c  fly  .  u 

(71(777)]  c  (z  +  ro> 

T  W 

2,2  2.  z  +  r 

c  +  (x  +  y  )  o 

/c2-(x2  +  y2)\  1 

(uy— vx) 

+  W 

4cfiy  v 

(c2+  (x^+  y  )]  c  (z  +  ro) 

c2  +  (x2  +  y2)  z  +  ro 

dw 
dT ' 


,2  2. 

(u  +  v  )  . 

z  +  r 


40  c 

2  ,  ,  2 

c  +  (x 


+  y*) 


(uy~  vx). 


(23) 


A  large  number  of  orthogonal  maps  are  such  that  the  longitude  and  latitude 
lines  on  the  earth  are  transformed  into  two  sets  of  mutually  perpendicular 
straight  lines.  In  these  x  is  a  function  of  X  alone  and  y  is  a  function  of  <j>  alone. 
Equivalently,  we  assume 


ax 

9y 


o. 


(24) 


8  X. 

(Actually  we  need  only  assume  the  first  *  0,  since  the  orthogonality  condition 
implies  *  0.  which  implies  0,  otherwise  *  0  andyy*  0  which  means 

y  is  not  a  function  of  X  or  0.  Conversely,  if  we  had  not  assumed  orthogonality, 

9  X 

Eqs.  (24)  trivially  imply  it  _  see  Eq.  [12].)  Assume  also  -yj-  >o  (the  case 
<  o  is  similar).  Then  Eqs.  (17)  become: 


Fx  ■  gj  +  (w  -  vej  tan  $)  (20cos  d  +  -jr). 


Fy  ■  ■g-j  +  ue  ^  tan  4  (2  ft  cos  $  +  -i)  +  , 
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2  2 

p  m  dw-u  +  v  —  2ft  cos  4  u. 

*  z  dt  r 

|*  §T"r  “u  (2°  C0B  *  +7>j' 

Here  e^  ■  sign  ("|y')"  *1‘ 

If  we  impose  the  further  restriction  upon  the  projection  that  it  be  conformal  we 
obtain  the  familiar  Mercator  Projection.  The  transformation  equations  in  this 
case  are: 

x  ■ 

y  -  faec*d  log  tan(-£+ ^  (26) 


or  the  inverse  set: 


\  *  x. 

$  ■  2  tan  *  e^—  w/2,  (27) 

r  ■  z  +  r  . 
o 

Equations  (25)  are  then  easily  written  entirely  in  the  new  system  by  substi¬ 
tuting 


tanh  y 

for 

sin 

sech  y 

for 

C08 

sinh  y 

for 

tan 

and  of  course 

z+ro 

for 

r. 

Another  large  class  of  maps  is  that  of  the  so-called  conic  projections.  These 
are  obtained  by  mapping  the  earth  onto  a  cone,  which  is  secant  to  the  earth 
(usually  such  that  the  cone's  axis  coincides  with  the  earth's  axis)  subject  to  the 
condition  that  points  on  the  same  latitude  line  map  into  points  at  the  same  height 
on  the  cone  and  points  on  the  same  longitude  line  map  into  points  on  the  same  ray 
of  the  cone.  The  cone  is  then  slit  along  a  ray  and  flattened  onto  a  plane. 
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The  transformation  can  be  expressed: 
p*  f  (4),  (f  an  arbitrary  monotonic  function). 


z  -  r-rQ. 

where  p  ,  a  are  polar  coordinates  on  the  x  -  y  plane,  that  is, 

f~2  2~ 

P *VX  +y  ' 

(*■  tan-1  {y  lx) . 

The  inverse  transformation  is  then 

4  ■  f-1  (p)  *  f_1  (v^T72) 

X.  ■  c  a  ■  c  tan-1  (y/x)  ,  (28) 

r  ■  z  +  r_. 

o  * 

Unfortunately,  such  a  map  is  not,  in  general,  orthogonal.  We  can  alter  our 
point  of  view  slightly,  however,  so  as  to  apply  the  equations  of  motion.  We  do  this 
by  thinking  in  terms  of  the  polar  coordinates  (  p,  a  )  instead  of  the  rectangular 
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coordinates  (x,  y),  so  that  "with  respect  to  p ,  a  "  the  transformation  is  orthogo¬ 
nal.  Indeed  we  are  now  in  the  situation  Eq.  (24)  with  x  and  y  replaced  by 
a  and  p,  respectively.  Thus,  if  we  now  look  at  radial  (e^  )  and  tangential  ( ea  ) 
directions  for  a  given  point,  P,  and  let  Fa  ,  F^  be  the  force  components  in 

A 
e/> 


X 

Figure  3.  Radial  and  Tangential  Directions 

these  directions  at  P  and  u,  v,  the  velocity  components,  then  Eqs.  (25)  apply, 
again  with  x  and  y  replaced  by  a  and  p  ,  respectively. 

Getting  back  to  rectangular  coordinates,  some  conic  projections  are  orthogo¬ 
nal  in  the  old  sense.  Applying  condition  Eq.  (12)  to  Eqs.  (26)  one  finds  the 
function  f  must  take  the  form 


P 


f  <tf)  -  k 


[tan(l+ f)] 


-1/c 


(29) 


In  this  special  case  the  transformation  also  turns  out  to  be  orthogonal  and  we 
have  in  fact  the  "Lambert"  confosmal  conic,  k  and  c  are  chosen  usually  to 
satisfy  specific  conditions,  such  as  true  scale  at  two  standard  parallels.  (Note 
that  the  stereographic  projection  is  an  instance  of  this  with  c  ■  1  and  k  ■  £ 
[that  is,  the  "c"  of  stereographic  projection]  since 

[•*n(i*D]'1TT&?- 


Of  course  this  was  to  be  expected  since  projection  on  a  plane  is  a  limiting  case 
of  projection  on  cones.) 

Thus,  Eqs.  (17)  apply  also  to  this  oft-used  class  of  projections,  and  take  a 
form  similar  to  that  for  stereographic  projection. 
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